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Preface 


This  report  contains  the  results  of  my  efforts  to 
develop  a  model  for  time-dependent  x-ray  transport  in  air. 
I  have  tried  to  present  a  complete  description  of  the 
derivation  of  the  equations,  the  assumptions  made,  and 
the  complexities  involved  in  trying  to  devise  a  set  of 
numerical  difference  equations  to  solve  the  model 
efficiently. 

There  are  several  people  to  whom  I  am  indebted  and 
owe  thanks.  First,  I  would  like  to  thank  my  advisor. 
Major  George  H.  Nickel,  for  without  his  help,  guidance, 
and  encouragement  I  would  not  have  completed  the  study. 

I  would  also  like  to  thank  Doctor  Charles  J.  Bridgman  and 
Captain  David  D.  Hardin  to  whom  I  frequently  went  for 
advice.  Lastly,  I  thank  my  wife,  Yuao-Ching,  for  her 
patience,  love,  and  understanding  throughout  the  period 
of  this  work. 

Robert  P.  Dickey 


(This  thesis  was  typed  by  Sharon  A.  Gabriel) 
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Abstract 


An  x-ray  transport  model  is  developed  which  gives  the 
streaming  radiation  flux,  thermal  radiation  flux,  and 
material  energy  of  air,  as  a  function  of  both  position  and 
time,  following  the  release  of  a  large  amount  of  x-ray 
energy  in  the  atmosphere.  The  model  is  intended  as  a 
compromise  between  the  simple  diffusion  theory  models  and 
the  accurate  transport  theory  models.  It  applies  to  the 
radiative  growth  which  takes  place  at  early  times  before 
hydrodynamic  motion  begins.  The  model  assumes  that  air  is 
heated  by  photoelectric  absorptions  to  a  burnout  temperature, 
the  kinetic  temperature  at  which  the  air  molecules  are 
completely  ionized.  The  opacity  of  air  is  assumed  very 
large  at  temperatures  below  the  burnout  temperature  and 
very  small  at  temperatures  above  the  burnout  temperature. 
Peierls'  integral  equation  is  used  to  calculate  the  flux 
of  thermal  radiation  which  builds  up  in  the  burned-out 
sphere  of  air.  The  integral  equation  considers  both  the 
nonlocal  and  retarded  time-dependent  characteristics  of 
the  thermal  radiation  field.  An  efficient  numerical  solution 
to  the  model  is  greatly  complicated  due  to  the  inverse  cubic 
temperature  dependence  of  the  opacity  and  the  retarded  time 
dependence  of  the  thermal  radiation  field.  A  numerical 
algorithm  is  presented  for  solving  the  model. 
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DEVELOPMENT  OF 


AN  X-RAY  TRANSPORT  MODEL 
BASED  ON  PEIERLS '  INTEGRAL  EQUATION 

I .  Introduct ion 

The  point  source  release  of  a  large  amount  of  x-ray 
energy  in  the  atmosphere  initiates  an  extremely  complex 
series  of  physical  processes  in  which  the  x-ray  energy 
liberated  by  the  source  heats  the  air,  and  the  air  then 
radiates  to  cool  itself.  The  competition  between  the 
various  modes  of  energy  transport  and  the  complexities  of 
the  radiation  properties  of  air  make  a  direct  analytic 
attack  infeasible  (Ref  1:2).  The  usual  approach  is  to  use 
transport  theory  and  obtain  numerical  solutions  to  the 
equation  of  transfer,  also  called  the  transport  equation. 

The  transport  equation  is  an  integro-differential 
equation  which  is  a  mathematical  statement  of  the  gains 
and  losses  of  photons.  Even  for  one-dimensional  systems 
with  spherical  symmetry,  the  transport  equation  expresses 
the  specific  intensity  of  radiation  as  a  function  of  four 
independent  variables:  position,  frequency,  direction,  and 
time.  A  typical  x-ray  transport  problem  might  involve  up 
to  50  or  more  radial  zones  and  perhaps  10  frequency  groups 
and  10  directions.  In  the  numerical  solution  to  this 
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problem,  each  combination  of  position,  frequency,  and 
photon  direction  requires  a  separate  mesh  point.  The 
calculation  is  entirely  nonlocal  in  that  photons  emitted 
at  any  one  mesh  point  at  a  specific  time  can  contribute  to 
the  number  of  photons  at  any  other  mesh  point  at  a  later 
time.  Clearly,  to  solve  this  problem  numerically  would 
require  extensive  computer  time  and  storage  capability. 

An  alternative  approach  to  solving  the  x-ray  transport 
problem  would  be  to  use  diffusion  theory.  Diffusion  theory 
assumes  that  the  radiative  flux  is  proportional  _o  the 
gradient  of  the  thermal  radiation  energy  density.  Although 
this  assumption  is  very  crude,  it  is  a  widely  used  approxi¬ 
mation  in  radiation  transport  calculations  because  of  its 
simplicity  (Ref  2:51).  Unlike  transport  theory,  diffusion 
theory  is  strictly  local  in  the  sense  that  the  number  of 
photons  at  any  one  mesh  point  can  only  be  influenced  by  the 
closest  neighboring  mesh  points.  If  the  photon  mean  free 
path  is  small  compared  to  the  dimensions  of  the  system, 
then  diffusion  theory  may  be  applicable.  However,  for 
most  optically  thin  spherical  systems  this  is  not  the  case 
and  diffusion  theory  does  not  apply. 

A  different  approach  to  the  x-ray  transport  problem 
would  be  to  solve  an  integral  form  of  the  transport  equation 
often  referred  to  as  "Peierls'  Equation"  (Ref  2:31).  Peierls' 
equation  solves  for  the  all-angle  thermal  radiation  flux  in 
a  system  which  scatters  photons  isotropically,  and  therefore 
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it  contains  two  less  independent  variables  which  describe 
direction  in  the  full  equation  of  transfer.  The  actual 
thermal  radiation  flux  may  or  may  not  be  isotropic.  A  form 
of  the  integral  equation- can  be  derived  for  time-dependent 
radiative  transfer  which  accounts  for  the  nonlocal  effect 
of  the  radiation  field.  An  efficient  numerical  solution  of 
Peierls*  equation  should  improve  upon  the  accuracy  and 
simplicity  of  diffusion  theory  models  without  using  as 
much  computer  time  and  storage  capability  as  the  full 
transport  theory  models. 

PURPOSE 

The  purpose  of  this  thesis  was  to  develop  a  model  that 
would  describe  the  radiative  transport  process  which  takes 
place  following  the  release  of  a  large  amount  of  x-ray 
energy  in  the  atmosphere.  The  intention  was  to  develop  a 
model  that  would  improve  upon  the  accuracy  of  diffusion 
theory  models  while  reducing  the  computer  time  needed  by 
full  transport  theory  models.  A  model  based  on  Peierls' 
integral  equation  will  be  developed  and  analyzed  to  see  if 
it  satisfies  the  above  criteria. 

SCOPE 

This  model  describes  one-dimensional,  t ime -dependent 
x-ray  transport,  in  spherical  geometry,  for  the  point  source 
release  of  a  large  amount  of  x-ray  energy  in  the  atmosphere. 
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The  model  applies  at  very  early  times  before  shock  formation 
and  hydrodynamic  expansion  takes  place.  Rather  than  treat 
the  photon  frequency  as  a  continuous  variable,  all  photons 
will  be  grouped  together  and  average  properties  assigned. 

ASSUMPTIONS 

The  following  assumptions  have  been  made  to  clarify 
and  simplify  the  x-ray  transport  model: 

1.  The  source  emits  x-rays  isotropically. 

2.  The  scattering  of  x-rays  is  isotropic. 

3.  The  density  of  air  is  uniform  at  its  ambient 
value  because  of  the  absence  of  hydrodynamic 
motion  during  early  time  periods. 

4.  The  opacity  of  air  has  an  inverse  cubic 
temperature  dependence. 

OVERVIEW 

This  report  contains  five  sections  and  one  appendix. 
Section  II  gives  a  theoretical  background  description  of 
the  physical  processes  taking  place  during  x-ray  transport. 
Based  on  the  underlying  physics,  the  equations  needed  to 
solve  the  model  are  developed  in  Section  III.  In  Section  IV, 
the  resulting  model  is  discussed  and  compared  to  an  existing 
diffusion  theory  model.  Finally,  the  conclusions  and 
recommendations  are  presented  in  Section  V.  Appendix  A 
describes  a  numerical  algorithm  which  could  be  used  to 
solve  the  model. 
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II.  Background 

This  section  presents  a  background  description  of 
what  physically  takes  place  when  a  large  amount  of  x-ray 
energy  is  released  in  the  air.  A  discussion  of  the  under¬ 
lying  physics  of  x-ray  transport  in  air  will  help  to  account 
for  assumptions  and  equations  presented  later  in  the 
development  of  the  model.  Much  of  the  discussion  comes 
from  the  work  of  Glasstone  (Ref  3),  and  Pomraning  (Ref  4), 
who  dealt  with  x-ray  transport  as  it  pertains  to  nuclear 
weapon  effects.  The  ionization  of  the  air  will  be  discussed 
first,  followed  by  an  estimation  of  the  time  required  to 
form  a  Maxwell-Boltzmann  distribution  of  electrons  and 
ions. 

IONIZATION  OF  AIR 

Initially,  an  intense  point  source  emits  x-rays  with 
several  KEVs  of  energy  into  the  atmosphere.  Close  to  the 
source,  the  x-rays  which  arrive  first  will  strip  off 
electrons  from  the  air  molecules  due  to  photoelectric 
absorptions.  As  more  electrons  are  stripped  from  the  air 
molecules,  the  air  begins  to  heat  up  until  it  reaches  what 
is  called  the  burnout  temperature.  The  burnout  temperature 
is  the  kinetic  temperature  ascribed  to  the  system  of 
completely  ionized  air  molecules  whose  velocity  distribution 
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is  Maxwellian.  Photoelectric  absorption  is  no  longer 
possible  in  the  completely  ionized  air  and  thus  an  upper 
limit  exists  on  the  amount  of  x-ray  energy  which  can  be 
deposited  in  a  given  volume  of  air.  The  air  opacity  is 
assumed  to  be  very  large  for  temperatures  below  the  burnout 
temperature  and  very  small  for  temperatures  above  the 
burnout  temperature  (Ref  4:220).  As  the  source  continues 
to  emit  x-rays,  the  spherical  region  of  completely  ionized 
air  around  the  source  expands  and  grows. 

Inside  the  burned-out  sphere,  the  heated  air  begins  to 
radiate  and  cool  down.  Electrons  that  were  stripped  from 
air  molecules  emit  bremsstrahlung  radiation.  Photons  and 
moving  electrons  interact  by  compton  and  inverse  coi..pton 
scattering.  As  the  kinetic  air  temperature  cools  below  the 
burnout  temperature,  the  opacity  of  air  increases,  and  more 
source  x-rays  are  then  absorbed  maintaining  the  air  at  the 
burnout  temperature.  The  net  result  of  these  various 
interactions  is  the  buildup  of  a  thermal  radiation  field  in 
the  burned-out  sphere.  For  systems  whose  dimensions  are 
on  the  order  of  the  mean  free  path  for  absorption  of  photons, 
the  actual  thermal  radiation  from  the  medium  will  be  of  the 
order  of  the  blackbody  value  (Ref  5:28).  Normally, 
"temperature"  implies  thermodynamic  equilibrium,  which 
means  both  kinetic  and  blackbody  radiation  equilibrium 
(Ref  5:28).  The  term  "kinetic  temperature"  has  been  used 
here  because  the  effective  thermal  radiation  temperature 
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will  initially  be  much  lower  than  the  kinetic  temperature 
of  the  air  particles,  resulting  in  a  two-temperature 
problem. 

Growth  of  the  burned-out  sphere  due  to  source  x-rays 
will  continue  until  the  source  strength  is  insufficient  to 
both  continue  ionizing  new  layers  of  air  and  maintain  the 
cooling  air  at  the  burnout  temperature.  This  is  the  end 
of  the  burnout  phase  and  further  growth  takes  place  through 
a  diffusion  process  (Ref  4:220).  At  this  time  the  ionized 
air  is  at  a  fairly  uniform  temperature,  although  temperature 
gradients  do  exist,  particularly  near  the  advancing 
radiation  front  (Ref  3:65).  The  thermal  radiation  field 
continues  building  up  in  the  burned-out  volume  of  air. 

The  absorption  of  the  thermal  radiation  in  the  surrounding 
cold  air  causes  the  diffusive  growth.  The  source  is  no 
longer  capable  of  keeping  the  air  hot,  and  as  the  air 
temperature  decreases,  the  photon  mean  free  path  decreases 
and  growth  of  the  ionized  sphere  becomes  less  rapid. 
Eventually,  the  effective  thermal  radiation  temperature 
and  the  material  kinetic  temperature  of  the  air  will 
approach  each  other  as  the  radiation  and  material  tend 
toward  equilibrium  (Ref  4:225). 

FORMATION  OF  A  MAXWELL- BOLTZMANN  DISTRIBUTION 

An  estimate  of  the  time  required  for  the  electrons 
and  ions  to  reach  a  Maxwellian  distribution  is  needed  to 
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justify  using  a  Maxwell-Boltzmann  equilibrium  expression 
for  the  material  energy  density  of  the  air.  Glasstone  and 
Lovberg  define  a  collision  as  a  coulomb  force  interaction 
which  results  in  a  deflection  of  more  than  90  degrees 
(Ref  5:83).  This  collision  can  be  a  single  interaction 
between  a  pair  of  charged  particles  (a  short-range  encounter) 
or  a  combined  net  effect  due  to  multiple  interactions 
between  a  single  particle  and  many  other  particles  (a 
long-range  encounter)  (Ref  5:83). 

To  determine  the  cross  section  for  a  short-range 
encounter,  Glasstone  defines  the  impact  parameter  for  a 
90  degree  scatter,  bQ  ,  as  the  distance  of  closest 
approach  between  two  particles  in  the  absence  of  any 
electrostatic  forces  (Figure  1).  From  Coulomb's  Law,  it 
is  found  that  the  particles  are  a  distance  2b^  apart  at 
the  point  of  closest  approach  for  a  90  degree  scatter,  and 
that  at  this  distance  the  mutual  coulomb  energy  is  equal  to 
the  center-of-mass  (or  relative)  kinetic  energy  of  the 
interacting  particles  (Ref  5:88).  For  the  case  of  two 
charged  particles  carrying  the  unit  electron  charge  e  , 

(Ref  5:88): 

W  -  or  b0  =  cm  (2.1) 
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Figure  1 


Short-Range  Coulomb  Interaction 
for  90  Degree  Deflection  (Ref  5:88) 
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where 


W 


relative  kinetic  energy  of  the  particles 
in  ergs 


e  =  electron  charge  in  statcoulombs . 


The  cross  section  for  the  short-range  encounter,  oc  ,  can 
be  estimated  as  the  area  of  a  disc  of  radius  bp  (Ref  5:88), 


o 
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irb  2 
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ire4 

4W2 


(2.2) 


In  analyzing  the  scattering  collision,  on  a  microscopic 
scale,  due  to  long-range  encounters,  every  charged  particle 
may  be  regarded  as  being  surrounded  by  an  atmosphere  having 
a  net  charge  of  the  opposite  sign  (Ref  5:84).  By  postulating 
that  the  charged  particles  have  a  Boltzmann  distribution, 
and  assuming  the  gradient  of  the  electric  potential  can  be 
expressed  by  Poisson's  equation  based  on  Coulomb's  Law,  the 
effective  radius  of  the  oppositely  charged  atmosphere 
surrounding  a  positive  ion  can  be  shown  to  be  (Ref  5:85): 
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kT 

irn  e  2 
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2.35  x  104 


T(Kev) 

ne(#/cm3) 


cm 


(2.3) 


where 

X^  *  Debye  shielding  (or  screening)  distance 

k  *  Boltzmann  constant 
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T  *  kinetic  temperature  of  both  electrons 

and  positive  ions  assuming  a  Maxwell - 

Boltzmann  distribution 

n  =  electron  number  density, 
e 


The  cross  section  for  a  long-range  encounter  is 
estimated  by  determining  the  momentum  change  that  results 
from  the  coulomb  force  at  the  distance  of  the  impact 
parameter  b  ,  operating  over  a  path  length  2b  ,  for 
a  time  2b/v  ,  where  v  is  the  velocity  of  the  scattered 
particle  (Ref  5:91).  The  particle  is  assumed  to  have  been 
scattered  through  a  large  angle,  about  90  degrees,  when  the 
average  of  the  square  of  the  total  momentum  change, 

(AP)2  ,  has  increased  to  the  point  that  it  is  roughly 
equal  to  the  square  of  the  initial  momentum  of  the 
scattered  particle,  P2  (Ref  5:92).  The  distance  traveled 
by  the  particle  is  equivalent  to  a  mean  free  path  for  the 
scattering  process,  and  it  is  possible  to  define  a  cross 
section  for  the  long-range  encounter,  ,  as  (Ref  5:93) 
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2ire4 
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where 

(Debye  shielding  distance) 

A  *  (Impact  parameter  for  a  90°  scatter, 

for  a  short-range  encounter) 


11 


W*  *  actual  energy  of  the  scattered  particle  taken 

to  be  3/2  kT  if  energy  distribution  is  Maxwellian. 

Substituting  Eq  (2.1)  and  (2.3)  for  bp  and  , 

respectively,  gives 


T2.35  x  1° ‘’(ItfuL 

h£l  J 


Assuming  the  ionized  atmosphere  consists  of  electrons  and 
nitrogen  ions,  with 


gives 


4.8  x  10~10  esu 

.4  kev  (the  burnout  temperature  of  air) 
5.76  x  1020  electrons/cm3 

1.9  x  10' 9  ergs 
9.6  x  10'10  ergs 


3.4x10"  cm2 


The  mean  collision  time  between  electron-electron 

scatters,  t  ,  can  now  be  estimated  as  (Ref  5:95) 

’  ee 


(2.5) 


12 


.  -.-.w  ..IHteu**'.  -«4  „ 


tee  '  infer  •  S'4  *  10"13  sec  (2-« 

e  d  e 

where  v  -  1.45  x  109  cm/sec  (average  electron  velocity). 

V 

Similarly,  the  mean  collision  time  between  ion-ion  scatters, 
t  ^  ,  can  be  estimated  as  (Ref  5:97) 


t .  . 
n 


n .  ojV . 
id  i 


6.0  x  10"  1 0  sec 


(2.7) 


where 

n^  Z  5.37  x  1019  ions/cm3 

v^  2  9.1  x  106  cm/sec  (average  ion  velocity). 


Finally,  the  mean  collision  time  between  electron-ion 
scatters,  tgi  ,  can  be  estimated  as  (Ref  5:97) 


ei 


n. 

l  d  e 


=  3.8  x  10 


1  2 


sec 


(2.8) 


The  time  required  for  several  collisions  between  electrons 
and  ions  is  so  short  compared  to  the  time  dependence  for 
source  emission  of  x-rays,  which  is  typically  on  the  order 
of  shakes  (10"  sec),  that  essentially  a  Maxwell - Boltzmann 
distribution  forms  instantaneously. 

An  interesting  comparison  can  be  made  by  estimating 
the  equilibration  time,  t  ,  for  thermal  photons  to  be 


13 


'P  c 


(2.9) 


where 


absorption  mean  free  path  of  air  at  burnout 

temperature 

speed  of  light 


For  X  :  50  meters  (Ref  4:225) 


t  =  1.67  x  10"  sec 

P 


While  the  air  particles  reach  a  Maxwell-Boltzmann  distribution 
on  the  order  of  10" 10  seconds,  the  thermal  photons  reach 
equilibrium  on  the  order  of  10'7  seconds,  giving  a  two- 
temperature  problem  with  the  thermal  radiation  field 
initially  not  in  equilibrium  with  the  particles. 
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III.  Development  of  X-Ray  Transport  Model 


This  section  derives  in  spherical  geometry  the  equations 
needed  to  describe  the  transport  of  a  large  amount  of  x-ray 
energy  released  in  the  air.  Three  questions  are  developed, 
one  each  for  the  streaming  radiation  flux,  thermal  radiation 
flux,  and  material  energy  of  air.  The  three  variables 
depend  not  only  on  time,  but  also  on  radial  position.  The 
equation  for  the  streaming  radiation  flux  describes  a  point 
source  that  has  been  exponentially  attenuated  and  spherically 
diverged.  A  form  of  Peierls'  integral  equation  will  be 
derived  to  calculate  the  thermal  radiation  flux,  and  the 
material  energy  of  air  will  be  determined  by  writing  a  rate 
equation  which  accounts  for  the  gains  and  losses  of  material 
energy. 

The  source  time  dependence  and  its  properties  will  be 
described  first,  followed  by  the  derivations  of  the  streaming 
radiation  flux,  thermal  radiation  flux,  and  material  energy 
of  air,  respectively.  The  end  of  this  section  will  show 
how  energy  is  conserved  in  the  model. 

SOURCE  TIME  DEPENDENCE 

Assume  a  point  source  of  radiation  of  source  strength 
Y  ,  and  a  time  dependence,  f  (t )  ,  given  by  Pomraning 

as  (Ref  6:3) 
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f(t) 


1  (ct+  g)ea  ~  t o ^ 
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(3.1) 


where 

a 

= 

parameter  for  source  growth 

3 

= 

parameter  for  source  decay 

to 

= 

parameter  for  the  time  when 

the  i 

maximum 

x-ray  energy  is  emitted  from 

the 

source 

Z 

= 

normalization  constant  such 

that 

the 

/  f(t)  dt  -  1  . 

0 

It  is  assumed  that  the  source  emits  a  Planck  blackbody 
spectrum  of  x-rays  into  an  infinite  sea  of  uniform-density 
air.  The  time  dependence  for  source  emission  of  x-rays 
has  the  following  properties  (Ref  6:3): 

(a)  Rises  like  eat  for  small  t  ; 

(b)  Falls  like  e"^t  for  large  t  ; 

(c)  Has  a  single  maximum  at  t  =  t  ; 

W)  fCt0)  -  \  . 


STREAMING  RADIATION  FLUX 

The  streaming  radiation  at  position  R  ,  time  t  , 
includes  only  those  virgin  x-rays  which  were  emitted  by  the 
source  at  an  earlier  time  and  arrived  untouched  at  position 
R  at  time  t  .  An  energy  balance  is  performed  for  the 
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streaming  radiation  density  in  3r  about  position  R 
and  in  at  about  time  t  (Ref  6:6). 


4ttR2Es  (R,t)  c  at 


/ 


-  4ttR2e  (R,t)  c  at 

S  i 


4wR23r  d[Es(R,t)] 


r+3r 


(3.2) 


Streaming  Radiation 

Streaming  Radiation 

Streaming  Radiation 

In 

Out 

Absorbed 

Here  c  is  the  speed  of  light  and  Es(R,t)  is  the  streaming 
radiation  density  at  position  R  ,  time  t  .  Eq  (3.2)  gives 
the  result  that  (Ref  6:7) 


£  [4TiR2Es(R,t)c]  =  4,tR2  [Es(R,t)] 


(3.3) 


The  streaming  radiation  emitted  by  the  source  at 
time  t  can  only  decrease  due  to  absorption  by  the  air 
as  it  streams  away  from  the  source.  Recall  the  opacity  of 
air  was  assumed  very  large  at  temperatures  below  the 
burnout  temperature  and  very  small  at  temperatures  above 
the  burnout  temperature.  It  is  now  assumed  that  the 
specific  temperature  dependence  of  the  opacity  is  given 
by  (Ref  6:15) 


<  (R,t) 


-  K 


B 

T(R,t) 
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(3.4) 
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1 


where 


k  (R ,  t ) 


K 


B 

T(R,t) 


the  opacity  of  air  at  position  R  , 
time  t  ,  which  depends  on  the  air 
temperature  at  position  R  ,  time  t  . 
the  opacity  of  air  at  the  burnout 
temperature  (a  constant), 
burnout  temperature  of  air 
the  kinetic  air  temperature  at 
position  R  ,  time  t 


The  time  rate  of  change  of  E  (R,t)  can  therefore  be 
expressed  as 


■k 


■ptc. 


tb 

Tim 


Es(R,t)c 


(3.5) 


where  p  is  the  density  of  air,  assumed  constant,  and 
<s  is  the  mean  opacity  of  air  at  the  burnout  temperature 
for  streaming  x-rays  being  emitted  by  a  Planckian  source. 

The  right  side  of  Eq  (3.5)  is  simply  the  rate  at  which 
streaming  radiation  density  is  absorbed  by  the  air. 

Equation  (3.3)  can  now  be  written  with  the  right  side 
of  Eq  (3.5)  substituted  in  for  the  time  rate  of  change  of 
the  streaming  radiation  energy  density,  giving 

TTKTtT  Es(R,t)c  (3.6) 


[4ttR2Es  (R,t)c]  =  -4ttR2p 
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Notice  that  Eqs  (3.7)  and  (3.9)  include  the  necessary 
time  retardation  which  accounts  for  the  finite  speed  of 
light  and  the  dependence  of  the  air  temperature  on  the 
path-length  variable  r  ,  in  the  exponential  attenuation 
term.  A  final  reduction  can  be  made  to  Eq  (3.9)  by 
recognizing  that 

Eg  (R,t  )c  =  Rs(R,t)  (3.10) 

where  Fs(R,t)  is  the  streaming  radiation  flux  at  position 
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R  ,  time  t  .  The  equation  for  the  streaming  radiation 
flux  can  now  be  written  in  final  form  as 


FsCR,t) 


Yf  (t-R/c) 
4ttR2 


B 


13 


T(r,t-'EI)J 


dr 


(3.11) 


The  streaming  radiation  flux  at  position  R  ,  at  time  t  , 
is  simply  the  source  strength  emitted  at  an  earlier  time, 
exponentially  attenuated  and  spherically  diverged  to 
arrive  at  position  R  ,  at  time  t 

THERMAL  RADIATION  FLUX 

The  thermal  radiation  flux  which  builds  up  in  the 
optically  thin,  heated  air  will  now  be  considered.  This 
consideration  must  include  both  the  nonlocal  and  retarded 
time  dependent  characteristics  of  the  thermal  radiation 
field.  The  development  of  the  equation  for  the  thermal 
radiation  flux  follows  an  analogous  development  by  Weinberg 
and  Wigner  for  an  integral  equation  for  the  transport  of 
monoenergetic  neutrons  (Ref  7:182). 

First,  a  differential  volume  element  in  the  shape  of 
a  flat  disc  is  positioned  at  R  ,  with  its  axis  in  the 
direction  from  R'  to  R  (Figure  2).  The  thermal 
radiation  being  emitted  at  position  R'  ,  at  time  t  , 
can  be  expressed  as 
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Figure  2.  Geometry  Associated  with  Thermal  Radiation  Flux 


where  "a"  is  the  radiation  constant  and  >ct  is  the 
mean  opacity  of  air  at  the  burnout  temperature  for  thermal 
x-rays.  In  order  to  arrive  at  position  R  ,  at  time  t  , 
photons  had  to  be  emitted  in  the  solid  angle  d (area)/4nS2 
at  an  earlier  time,  t  -  S/c  ,  where  S  is  the  path 
length  from  R'  to  R  ,  and  d(area)  is  the  differential 
area  of  the  base  of  the  disc  (Figure  2).  The  probability 
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that  photons  arrive  at  the  disc  without  undergoing  further 


attenuation  is 
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(3.13) 


Since  the  height  of  the  disc  is  equal  to  dh  ,  the  sum  of 
the  path  lengths  of  the  photons  in  the  disc  at  time  t  , 
due  to  photons  emitted  earlier  at  position  R'  is 


T(tr,t-S/c) 


caT(R',t-S/c)4  tett 
4itS2 


S 

"/  PKt 


T(k',t-S'/c) 
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(3.14) 


Equation  (3.14),  divided  by  the  volume  of  the  disc, 
d(area)dh  ,  is  the  contribution  to  the  thermal  radiation 
flux  at  position  R  ,  at  time  t  ,  due  to  radiation  emitted 
earlier  from  position  R'  .  The  total,  all-angle  thermal 
radiation  flux,  at  position  R  ,  at  time  t  ,  Ft(R,t)  , 
will  be  the  integral  of  Eq  (3.14)  over  the  entire  volume, 
namely 


Ft(R,t)  =  /dV OL  PKt 


caT(R',t-S/c) 
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(3.15) 
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Equation  (3.15),  integral  in  character,  describes 
time  dependent  radiative  transfer  in  a  system  which  isotropi¬ 
cally  scatters  photons.  The  actual  flux  itself  will  not  be 
isotropic,  especially  near  the  boundaries.  The  equation  is 
of  a  form  commonly  referred  to  as  Peierls'  equation  (Ref  2:31). 

The  fact  that  the  opacity  depends  on  one  of  the  dependent 
variables,  the  kinetic  air  temperature,  greatly  complicates 
the  solution  to  Eq  (3.15).  It  will  be  useful  to  define  an 
effective  temperature  for  the  thermal  radiation  field, 

T^^  ,  from  the  relationship 

oTeff(R,t)4  =  Ft(R,t)  (3.16) 

where  G  is  the  Stefan -Boltzmann  constant. 

MATERIAL  ENERGY  DENSITY 

The  material  energy  density  of  air  is  simply  the 
kinetic  energy  of  the  ions  and  electrons.  In  a  given 
molecular  system,  the  exchange  of  energy  between  colliding 
ions  and  electrons  causes  the  presence  of  a  Maxwell- 
Boltzmann  distribution  of  energies.  Earlier,  it  was  shown 
that  a  Maxwell -Boltzmann  distribution  of  electrons  and  ions 
forms  almost  instantaneously  following  the  release  of  a 
large  amount  of  x-ray  energy  in  the  air.  Therefore,  based 
on  the  assumption  that  the  molecular  air  is  in  equilibrium 
with  a  Maxwell -Boltzmann  distribution  of  energies,  the 
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-t*: 


where  Em(R,t)  is  the  material  energy  density  of  air,  and 
Cv  is  the  heat  capacity  of  air,  assumed  constant  (Ref  4:223). 
The  real  unknown  is  the  kinetic  air  temperature,  since  the 
material  energy  density  is  known  if  the  temperature  is  known. 

Three  processes  are  taking  place  which  affect  the 
temperature  of  air.  The  temperature  of  air  is  increased 
when  either  streaming  radiation  or  thermal  radiation  is 
absorbed,  and  is  decreased  when  the  air  radiates  to  cool 
itself.  An  equation  is  written  for  the  time  rate  of  change 
of  the  material  energy  density  due  to  these  gains  and  losses. 


(3.18) 

An  iterative  process  will  be  required  to  numerically  solve 
for  the  unknown  air  temperature  in  Eq  (3.18). 
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CONSERVATION  OF  ENERGY 


This  x-ray  transport  model  was  developed  from  the  basic 
principle  of  energy  conservation.  At  any  time  t  ,  the 
total  streaming  radiation  energy,  thermal  radiation  energy, 
and  material  energy,  contained  in  the  volume  of  air,  must 
equal  the  total  energy  emitted  by  the  source  up  to  time  t 
When  streaming  radiation  or  thermal  radiation  is  attenuated 
by  absorption,  the  material  energy  of  air  must  increase. 
Likewise,  when  the  air  radiates  to  cool  itself,  thermal 
radiation  energy  must  increase  while  material  energy  must 
decrease. 


On  inspecting  the  equations  for  the  streaming  radiation 
flux,  thermal  radiation  flux,  and  material  energy  density 
of  air,  it  is  not  obvious  that  energy  is  conserved.  In 
order  to  show  that  energy  is  conserved,  the  equations  for 
the  streaming  radiation  and  the  thermal  radiation  are 
expressed  as  energy  densities  instead  of  fluxes.  The  three 


energy  density  equations  are  now 
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where  the  time  rate  of  change  of  the  material  energy  density 
was  replaced  by  a  backward  difference  equation  and  Et(R,t)  , 
which  is  the  thermal  radiation  density  at  position  R  and 
time,  defined  from  the  relationship 


|Et(R,t)  =  Ft(R,t)  (3.22) 

Conservation  of  energy  is  now  determined  by  integrating 
the  three  equations  over  the  entire  volume  and  adding  them 
up.  Integrating  and  adding  up  the  left-hand  sides  of 
Eqs  (3.19),  (3.20),  and  (3.21)  gives  the  total  streaming 
radiation  energy,  thermal  radiation  energy,  and  material 
energy  in  the  volume  at  time  t  .  On  integrating  and  adding 
up  the  right-hand  sides  of  Eqs  (3.19),  (3.20),  and  (3.21), 
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it  is  found  that  the  exponential  attenuation  portion  of 
Eq  (3.19),  integrated  over  the  volume,  is  cancelled  by  the 
material  energy  gains  due  to  streaming  radiation  absorption, 
terms  ©  and  ©  Eq  (3.21).  Similarly,  the  emission  and 
exponential  attenuation  portions  from  the  thermal  radiation 
energy,  term  ( 1 )  Eq  (3.20),  are  cancelled  by  terms  ©  , 

©  ,  and  ©  ,  in  Eq  (3.21)  when  integrated  over  the  entire 
volume.  What  remains  after  adding  up  all  terms  is  only 
the  source  strength  emitted  up  to  time  t  in  term  ©  of 
Eq  (3.19).  Energy  is  conserved  since  the  total  streaming 
radiation  energy,  thermal  radiation  energy,  and  material 
energy  is  equal  to  the  source  strength  emitted  up  to 
time  t 
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IV.  Results  and  Discussion 


In  this  section  the  x-ray  transport  model  developed 
from  Peierls'  integral  equation  will  be  analyzed.  Recall, 
the  intention  was  to  develop  a  model  that  would  improve 
upon  the  accuracy  of  diffusion  theory  models  while  reducing 
the  computer  time  required  of  full  transport  theory  models. 
Based  on  that  criterion,  the  model  developed  in  this  study 
will  be  compared  to  an  existing  diffusion  theory  model  by 
Pomraning  (Ref  6).  A  qualitative  measure  of  the  accuracy 
will  be  made  by  comparing  the  basic  theory  and  assumptions 
on  which  each  model  is  based.  The  model  developed  in  this 
study  will  then  be  numerically  analyzed  to  see  if  the 
model  lends  itself  to  an  efficient  numerical  solution. 

COMPARISON  TO  POMRANING' S  DIFFUSION  THEORY  MODEL 

Pomraning  develops  a  relatively  straightforward  model 
for  the  release  of  a  large  amount  of  x-ray  energy  in  the 
atmosphere  by  assuming  two  distinct  phases  for  the  growth 
of  the  ionized  sphere  of  air  (Ref  4:220).  The  initial 
phase,  or  burnout  phase,  consists  of  the  complete  ioniza¬ 
tion  and  burnout  of  the  air,  followed  by  a  diffusive  growth 
phase  which  begins  when  the  source  of  x-rays  can  no  longer 
maintain  the  air  at  the  burnout  temperature.  The  final 
conditions  that  exist  at  the  end  of  the  burnout  phase  become 


28 


”5* 


the  initial  conditions  for  the  start  of  the  diffusive 
growth  phase.  In  reality  there  are  not  two  distinct  phases, 
but  a  smooth  blending  from  one  phase  to  the  next.  The 
model  developed  in  this  study  does  not  separate  the  two 
phases,  but  allows  them  to  blend  smoothly  together. 

In  Pomraning’s  burnout  phase,  he  assumes  the  sphere  to 
be  optically  thin  and  therefore  he  ignores  spatial  gradients 
and  writes  rate  equations  for  the  volume  as  a  whole  (Ref  4:221). 
He  advances  a  radiation  front  with  an  associated  burnout 
radius.  Inside,  the  burnout  radius  is  completely  ionized  air 
at  the  burnout  temperature,  while  outside  the  burnout  radius 
is  ambient  air.  The  model  developed  in  this  study  also 
assumes  an  optically  thin  sphere,  but  spatial  gradients  are 
not  ignored.  The  thermal  radiation  flux,  streaming  radiation 
flux,  and  kinetic  air  temperature  are  all  determined  as  a 
function  of  both  position  and  time.  According  to  Glasstone 
(Ref  3:65),  temperature  gradients  do  exist,  particularly 
near  the  advancing  radiation  front.  If  temperature  gradients 
exist,  then  thermal  radiation  gradients  would  also  exist,  and 
a  numerical  solution  to  this  model  would  give  a  good  indication 
of  the  magnitude  of  these  gradients.  Since  the  kinetic  air 
temperature  is  determined  as  a  function  of  position  and  time, 
it  is  no  longer  necessary  to  advance  a  radiation  front  with 
an  associated  burnout  radius.  The  spatial  dependence  of 
the  air  temperature  will  indicate  the  volume  of  air  that  is 
completely  ionized,  as  well  as  the  regions  that  are  only 
partially  ionized. 


29 


After  completion  of  the  burnout  phase,  Pomraning  uses 
diffusion  theory  to  model  further  growth  of  the  ionized 
sphere  of  air.  To  obtain  diffusion  theory,  he  assumes  the 
radiative  flux  is  proportional  to  the  gradient  of  the 
thermal  radiation  energy  density  (Ref  4:223).  Diffusion 
theory  is  most  accurate  when  the  dimensions  of  the  system 
are  large  compared  to  the  photon  mean  free  path.  For  a 
spherical  system  of  ionized  air  with  a  radius  in  the  tens 
of  meters,  all  points  in  the  system  are  within  a  few  mean 
free  paths.  Pomraning  also  assumes  that  the  kinetic  air 
temperature,  the  thermal  radiation  density,  and  the  streaming 
radiation  density  are  all  space  independent  within  the  sphere 
of  ionized  air,  and  zero  elsewhere  (Ref  4:224).  Instead  of 
the  diffusion  equation,  this  model  was  developed  around 
Peirels'  integral  equation.  Peirels'  equation  gives  spatial 
dependence  while  accounting  for  the  time  retardation  and 
nonlocal  characteristics  of  the  radiation  field.  The 
accuracy  of  Peierls'  equation  is  not  dependent  on  the  photon 
mean  free  path  or  the  nearness  of  a  calculation  to  a 
boundary. 

NUMERICAL  ANALYSIS 

Before  proceeding  directly  into  the  numerical  analysis, 
some  introductory  comments  will  be  made  on  the  numerical 
approach  to  solving  this  model.  The  initial  idea  for  using 
an  integral  equation  of  Peierls'  form  came  from  the 
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similarity  of  Peierls*  integral  equation  to  the  integral 
from  of  the  neutron  transport  equation.  Eugene  Skluzacek, 
in  his  doctoral  dissertation  (Ref  8),  presented  a  very 
efficient  numerical  scheme  for  solving  the  integral  form 
of  the  neutron  transport  equation.  His  efficient  numerical 
solution  relied  on  an  "interpolation  coefficient"  developed 
by  his  advisor,  George  K.  Nickel.  This  "interpolation 
coefficient"  (Ref  8:53)  was  used  in  the  integration  over 
S  to  express  the  flux  as  a  function  of  S  rather  than 
R  (see  Figure  2,  page  21).  The  initial  plan  was  to  apply 
this  same  numerical  scheme  to  the  integration  over  S  in 
Peierls'  integral  equation  for  the  thermal  radiation  flux. 

Two  very  important  differences  exist  between  Peierls' 
integral  equation  and  the  integral  form  of  the  neutron 
transport  equation.  These  two  differences  greatly  complicate 
using  the  same  numerical  scheme.  The  first  difference  is 
in  the  cross  sections.  The  neutron  cross  section  is 
uniform  between  radial  mesh  boundaries,  while  the  photon 
opacity  is  non-uniform.  Not  only  is  the  photon  opacity 
non-uniform,  but  it  also  depends  on  one  of  the  dependent 
variables,  namely  the  kinetic  air  temperature.  The 
second  difference  comes  in  the  time  dependence  of  the 
neutron  transport  problem  compared  to  the  x-ray  transport 
problem.  In  the  case  of  neutron  transport,  the  time 
dependent  problem  can  be  transformed  into  a  stationary 
problem  by  assuming  an  exponential  growth  rate  and  solving 
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the  resulting  eigenvalue  problem.  For  the  x-ray  transport 
model  developed  in  this  study,  the  time  dependence  -- 
more  specifically,  the  retarded  time  dependence  --  must  be 
numerically  dealt  with  to  obtain  an  accurate  solution  to 
the  model. 

The  equation  for  the  thermal  radiation  flux. 


r  B 

VR>t)  =  dv0L  pKt|-T"(R7',t - 


57c 7 


caT  (R',t-S/c) 


4irS2 


13 


T(R'vt"570 


dS' 


(4.1) 


is  now  re-examined  to  determine  the  effect  of  the  retarded 
time  dependence,  and  the  opacity  temperature  dependence,  on 
a  numerical  solution.  To  avoid  a  discontinuity  which 
arises  when  R'  =  R  ,  S  is  chosen  as  the  variable  of 
integration.  The  differential  volume  element  becomes 


dVOL  =  S2sin  e  de  d$  dS 


(4.2) 


where  e  is  the  angle  between  R  and  S  .  Since  the 
problem  is  symmetric  in  the  azimuthal  angle,  the  $ 


32 


9 


rm 


integration  can  be  evaluated  to  give 

dVOL  =  2uS2sin  e  de  dS  (4.3) 

The  integration  over  6  can  easily  be  handled  if  a  trans¬ 
formation  is  made  to  an  integration  over  p  =  cos  e  , 
since 


/  sin  e  def(e)  =  /  diif(y) 

0  -1 

Integrals  of  this  type  are  accurately  approximated  by 
Gaussian  quadrature,  giving 


(4.4) 


+1  L 

/  dyf  (u )  2  l  «  f(p.)  (4.5) 

-1  l-l 

where 

=  quadrature  weights 

u  =  direction  cosines. 
i 

The  volume  integral  has  now  been  reduced  by  azimuthal 
symmetry  and  Gaussian  quadrature  to  a  weighted  sum  of  L 
one-dimensional  integrals  over  S  ,  giving 
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In  the  neutron  transport  model,  the  exponential  growth  can 
be  incorporated  into  the  exponential  attenuation  term  by 
defining  a  new  cross  section  (Ref  8:17).  Since  the  neutron 
cross  section  is  uniform  between  radial  shells,  it  is 
possible  to  incorporate  the  "interpolation  coefficient"  for 
the  flux,  and  then,  by  a  change  of  variables,  get  an 
integral  of  the  form 

b  n  - 

/  Zn  e  dZ  (4.7) 

a 

This  integral  is  just  the  partial  gamma  function  which  has 
an  analytic  solution. 

In  Equation  (4.6)  for  the  thermal  radiation  flux,  the 
opacity  depends  on  the  temperature.  Since  the  temperature 
depends  on  position  and  time,  the  exponential  integral 
cannot  be  solved  analytically  in  its  present  form.  If  the 
kinetic  air  temperature  was  defined  on  the  shell  boundaries 
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and  varied  linearly  with  R  between  shell  boundaries, 
then  the  integration  over  S  would  require  interpolation 
in  both  position  and  time  for  the  air  temperature.  The 
crossing  of  shell  boundaries  at  retarded  times  during  an 
integration  over  S  is  not  likely  to  coincide  to  times 
when  the  kinetic  air  temperature  has  previously  been 
determined.  The  ’’interpolation  coefficients”  require 
modifications  to  account  for  the  dependence  of  the  air 
temperature,  and  hence  the  opacity,  on  position  and  time. 

The  original  difference  scheme  which  worked  very  efficiently 
for  the  transport  of  neutrons  was  found  to  be  too  difficult 
to  apply  efficiently  to  solving  this  x-ray  transport 
equation.  For  this  reason,  a  numerical  solution  to  this 
model  was  not  reached.  The  equations  defining  the  streaming 
radiation  flux  and  material  temperature  of  air  both  lend 
themselves  to  more  efficient  numerical  solutions  and  do 
not  introduce  any  specific  numerical  complications. 
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r.  Conclusions  and  Recommendations 

A  theoretical  x-ray  transport  model  has  been  developed 
to  describe  the  radiative  transport  process  which  takes 
place  following  the  release  of  a  large  amount  of  x-ray 
energy  in  the  atmosphere.  The  model,  based  on  a  form  of 
Peierls'  integral  equation,  was  intended  to  be  a  compromise 
between  simple  diffusion  theory  models  and  accurate  transport 
theory  models.  The  conclusions,  drawn  from  the  analysis  in 
Section  IV,  will  be  presented  first,  followed  by  the 
recommendations . 

CONCLUSIONS 

Radiative  transport,  as  described  with  a  form  of 
Peierls'  integral  equation,  accurately  models  the  physical 
processes  which  take  place  when  a  large  amount  of  x-ray 
energy  is  released  in  the  atmosphere.  The  integral  approach 
to  solving  for  the  thermal  radiation  flux  accurately  takes 
into  account  the  nonlocal  and  retarded  time  dependent 
characteristics  of  the  radiation  field.  Since  diffusion 
theory  is  limited  in  its  ability  to  treat  the  streaming 
radiation  flux,  two  distinct  phases  are  required  to  model 
the  radiative  transfer  based  on  diffusion  theory.  The 
model  developed  in  this  study  provides  a  smooth  transition 
from  growth  of  the  ionized  sphere  due  to  streaming  radiation 


to  growth  due  to  thermal  radiation.  Unlike  diffusion 
theory  models,  the  accuracy  of  a  solution  to  this  model  is 
not  affected  by  the  nearness  of  a  boundary  or  the  size  of 
the  system  in  comparison  to  the  photon  mean  free  path. 
Instead  of  ignoring  spatial  gradients,  this  model  calculates 
the  spatial  gradients,  giving  more  detail  and  resolution 
to  the  resulting  values  of  the  streaming  radiation  flux, 
thermal  radiation  flux,  and  material  temperature  o."  air. 

From  a  theoretical  standpoint,  this  model  should  be 
considerably  more  accurate  than  diffusion  theory  models. 

A  numerical  solution  to  this  x-ray  transport  model  is 
greatly  complicated  due  to  time  retardation  and  the 
temperature  dependence  of  the  opacity.  A  numerical 
difference  scheme  that  had  worked  very  efficiently  for  the 
integral  form  of  the  neutron  transport  equation  was  found 
to  require  too  many  approximations  to  account  for  the 
opacities  temperature  dependence  and  time  retardation. 

These  approximations  served  to  increase  the  complexity  of 
the  numerical  solution,  while  decreasing  its  efficiency. 

The  final  conclusion  was  that  an  efficient  numerical 
solution  to  this  model  could  not  be  obtained  with  the 
originally  planned  numerical  difference  scheme.  There 
was  not  enough  time  to  devise  an  alternative  difference 
scheme  that  would  solve  the  model  efficiently,  without 
making  simplifying  assumptions  and  reducing  the  model’s 
accuracy.  Given  more  time,  it  is  possible  that  an 
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efficient  numerical  solution  to  the  model  might  be  worked 
out . 

RECOMMENDATIONS 


The  following  recommendations  present  the  author's 
view  concerning  extensions  of  this  research  and  possible 
applications  for  this  x-ray  transport  model. 

1.  Appendix  A  presents  an  algorithm  for  solving  a 
simplified  version  of  the  model,  where  the  kinetic  air 
temperature  at  any  time  t  ,  is  defined  between  shell 
boundaries  and  assumed  uniform  between  shell  boundaries. 
Continued  work  is  needed  to  devise  an  algorithm  and  a  set 
of  difference  equations  for  the  more  difficult  problem 
where  the  kinetic  air  temperature  is  defined  on  the  shell 
boundaries  and  assumed  to  vary  linearly  between  shell 
boundaries.  From  the  set  of  difference  equations,  a 
numerical  solution  can  be  obtained  and  a  measure  of  the 
model's  accuracy  and  efficiency  determined  by  comparison 
to  other  x-ray  transport  benchmark  calculations. 

2.  Throughout  this  model,  it  was  assumed  that  the 
opacity  of  air  had  an  inverse  cu^ic  temperature  dependence. 
This  assumption,  taken  from  Pomraning's  model  (Ref  4:223), 
turned  out  to  be  a  critical  assumption,  as  it  caused  most 
of  the  numerical  complications.  It  is  recognized  that  the 
photoelectric  cross  section  is  inversely  proportional  to 
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the  cube  of  the  incoming  photon  energy,  but  the  fact  that 
the  opacity  has  an  inverse  temperature  dependence  is  not 
obvious.  Although  there  is  no  reason  to  doubt  Pomraning’s 
assumption,  the  question'  has  been  raised  as  to  justification 
of  the  assumption.  Continued  work  on  the  numerical  solution 
to  this  model  should  include  justifying  this  assumption 
and/or  investigating  other  alternative  approaches  to 
determining  the  opacity  of  air. 

3.  If  an  efficient  numerical  solution  can  be 
developed  for  this  model,  then  consideration  should  be 
given  to  using  the  output  from  this  model  as  the  input 
to  a  hydrodynamic  model  which  develops  and  propagates 
the  shock  wave. 
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Appendix  A 


Numerical  Algorithm  for  Solving  the  Model 


This  appendix  presents  an  algorithm  for  a  numerical 
solution  to  this  x-ray  transport  model.  The  algorithm  is 
not  intended  to  be  a  step-by-step  computer  implementation 
of  a  set  of  difference  equations,  but  merely  serves  as  an 
outline  to  indicate  the  logic  and  direction  of  flow  for  a 
computer  program  to  follow.  The  algorithm  will  be  preceded 
by  a  general  description  of  a  typical  radial  mesh,  to 
include  the  radial  dependence  of  the  various  quantities 
being  solved  for. 


TYPICAL  RADIAL  MESH 

A  typical  radial  mesh  would  consist  of  a  series  of 
concentric,  spherical  shells  whose  boundaries  can  be 
defined  by  the  relationship  (Figure  A-l), 


ri.i  ■  ri  *  4r 


(A.l) 


Both  the  streaming  radiation  flux  and  the  thermal  radiation 
flux  would  be  defined  on  the  shell  boundaries.  It  is 
assumed  that  between  specific  radii  the  streaming  and  thermal 
radiation  fluxes  vary  linearly  with  r  .  The  algorithm  is 
much  more  difficult,  and  has  not  been  worked  out  for  the 
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case  when  the  kinetic  air  temperature  is  defined  on  the 
shell  boundaries,  and  is  assumed  linear  with  r  between 
shell  boundaries.  This  algorithm  will  treat  the  simpler 
case  of  the  kinetic  air  temperature  being  determined  between 
shell  boundaries,  and  assumed  uniform  between  shell  boundaries. 
Figure  A-l  depicts  a  radial  mesh  along  with  several  of  the 
Gaussian  quadrature  rays  that  will  be  needed  for  the  volume 
integration  to  determine  the  thermal  radiation  flux. 

A  typical  x-ray  transport  problem  might  require  50 
radial  shells,  each  shell  being  1.5  meters  apart.  To  trace 
the  radiative  growth  from  time  t=0  to  t  =  1.0  x  10“6  sec 
would  require  at  least  200  time  steps.  Usually,  the  time 
steps  are  no  longer  than  the  time  it  takes  a  photon  to  travel 
the  width  of  a  radial  shell.  The  volume  integral  for  the 
thermal  radiation  flux  would  require  approximately  ten 
(L  *  10)  Gaussian  quadrature  rays  to  get  acceptable  convergence 
of  the  0  integration.  For  50  radial  shells  and  200  time 
steps,  a  minimum  of  30,000  memory  locations  would  be  required 
for  the  three  dependent  variables. 

NUMERICAL  ALGORITHM 

The  following  algorithm  lists  in  order  the  steps  to  be 
taken  in  solving  this  x-ray  transport  model.  The  algorithm 
is  for  the  case  when  the  time  increments,  At  ,  are  exactly 
equal  to  Ar/c 
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STEP  1 

STEP  2 

STEP  3 


STEP  4 

STEP  5 

STEP  6 


Read  in  all  the  applicable  constants  and  parameters. 
This  includes  Y  ,  a  ,  B  ,  t„  ,  c  ,  k  ,  TD  ,  k_ 

O  X5  ^ 

p  ,  o  ,  a  ,  and  . 

Normalize  the  time  dependence  and  determine  a 
value  for  Z  based  on  a  ,  B  ,  and  t  .  An 
efficient  two-point  quadrature  rule  is  presented 
by  Pomraning  for  this  normalization  (Ref  6:23). 

Read  in  the  initial  conditions.  The  air  temperature, 
for  all  mesh  regions  and  times  can  be  initially  set 
to  standard  atmospheric  temperature.  The  thermal 
radiation  flux,  for  all  radii  and  times,  can 
initially  be  set  to  zero.  The  streaming  radiation 
flux  can  be  set  equal  to  zero  for  all  radii  (except 
r  =  0)  and  times.  Select  a  value  for  Ar  ,  At  , 
and  calculate  the  mesh  radii  and  times  from 


rl  "  0 


ri+l 


h  =  0 


■n+l 


Increment  the  time  index 

Set  the  radial  index,  i 
index,  n 


Calculate  F  n  from  Eq 
si 


r-_  +  Ar 

t  +  At 
n 

,  n  . 

,  equal  to  the  time 


(3.11). 


(A.  2 
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STEP  7 


Calculate  T?  from  Eq  (3.18). 

STEP  8  Calculate  F  n  from  Eq  (3.15). 

zi 

STEP  9  Repeat  STEPS  6,  7,  and  8,  keeping  the  time  index, 

n  ,  constant  while  reducing  the  radial  index,  i  , 
by  one  each  time  until  Fg  ,  T  ,  and  F  have 
been  calculated  for  all  radii  i  ,  less  than  or 
equal  to  n  .  Then,  go  back  to  STEP  4  and  increment 
the  time  index,  and  start  the  set  of  calculations 
over.  The  calculation  continues  until  the  final 
value  of  the  time  is  reached. 

ADDITIONAL  COMMENTS 

At  late  times,  the  calculation  for  the  streaming 
radiation  flux  can  be  omitted  when  the  source  strength, 

Yf(t)  ,  is  negligibly  small. 

The  time  rate  of  change  of  the  air  temperature  in 
Eq  (3.18)  can  be  replaced  by  a  backward  difference  equation. 

The  unknown  air  temperature  will  have  to  be  solved  with  an 
iterative  procedure.  One  approach  is  to  use  a  Bolzano  search 
which  halves  the  interval  of  uncertainty  each  iteration. 

While  integrating  on  one  of  the  Gaussian  quadrature 
rays  from  S  =  0  to  S  ,  in  the  calculation  of  the  thermal 
radiation  flux,  the  air  temperature  will  change  when  either 
a  new  mesh  region  is  entered  or  a  retarded  time  is  reached 
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for  which  a  different  air  temperature  applies.  A  test 
should  be  made  which  terminates  the  integration  when  either 
the  contributions  are  negligibly  small  or  a  retarded  time  is 
reached  that  is  negative. 
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